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DIRECT  RUNOFF  MVDROGRAPH 
PARAFETERS  VERSUS  URBANIZATION 

By  David  L.  Gundlach,^  A.M.  ASCE 


Introduction 

Various  rainfall-runoff  models  are  based  on  the  development  of  unit  hydro¬ 
graphs,  loss  rate  functions,  and  routing  criteria.  With  models  of  this  type, 
characteristics  used  to  define  the  unit  hydrograph,  loss  rate,  and  routing  criteria 
need  to  be  modified  to  predict  runoff  that  would  occur  because  of  future 
development  within  a  watershed.  Certain  aspects  of  this  problem,  particularly 
changes  in  peak  flow  and  lag  time  due  to  urbanization,  have  been  treated  previously 
(1,2,4,5,6).  It  is  the  aim  of  this  note  to  present  additional  information  regarding 
the  modification  of  unit  hydrograph  characteristics  due  to  increased  urbanization 
and  to  introduce  techniques  which  can  be  utilized  in  a  practical  solution. 

Effects  of  Urbanization 


A  multiple  regression  analysis  based  on  15  flood  hydrograph  reconstitutions 
in  the  vicinity  of  Philadelphia,  utilized  in  the  preliminary  report,  “Metropolitan 
Chester  Creek  Basin,  Pennsylvania,”  Department  of  the  Army.  Philadelphia 
District,  Corps  of  Engineers,  January  1976,  resulted  in  the  following  expressions 
(see  Table  1): 


(TC  +  R)  =  19.46  f 


( TC )  =  12.98  / 


/  DA  \° : 


Eqs.  1  and  2  relate  the  direct  runoff  hydrograph  parameters,  TC  and  R  to 
physiographic  characteristics  of  the  drainage  basin,  in  which  TC  =  the  time 
from  the  end  of  effective  rainfall  to  the  inflection  point  on  the  recession  limb 
of  the  hydrograph,  in  hours;  R  =  the  ratio  of  the  discharge  at  the  inflection 
point  on  the  recession  limb  of  the  hydrograph  to  the  rate  of  change  of  discharge 
at  that  point,  in  hours;  /  =  the  percentage  of  impervious  surface  within  a  watershed 

Note. — Discussion  open  until  February  I,  1977.  To  extend  the  closing  date  one  month, 
a  written  request  must  be  filed  with  the  Editor  of  Technical  Publications,  ASCK.  This 
paper  is  part  of  the  copyrighted  Journal  of  the  Irrigation  and  Drainage  Division.  Proceedings 
of  the  American  Society  of  Civil  Engineers,  Vol.  102,  No.  IR3,  September.  1976.  Manuscript 
was  submitted  for  review  for  possible  publication  on  October  23,  1975. 

'  Hydr.  Engr.,  Hydrologic  Engrg.  Center,  Corps  of  Engrs.,  U.S.  Dept,  of  the  Army. 
Davis,  Calif. 
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(a  measure  of  urbanization);  DA  =  drainage  area,  in  square  miles;  and  5  = 
the  average  channel  slope  between  the  points  I (T/t  and  90%  of  the  distance 
upstream  from  the  gage  or  outflow  point  to  the  watershed  boundary,  in  feet 
per  mile.  Although  the  correlation  was  improved  in  subsequent  work  when 
specific  physiographic  and  meteorological  characteristics  were  combined  (4). 
the  information  outlined  in  Table  I  is  sufficient  to  illustrate  the  following 
techniques. 

TABLE  1.— I  Results  of  Multiple  Regression  Analysis  In  Which  Direct  Runoff  Hydrograph 
Characteristics,  TC  +  R  and  TC,  are  Related  to  Physiographic  Characteristics  of 

Drainage  Basin 


Coefficient 

Standard 

Correlation 

of 

Equation 

error  of 

coefficient. 

determination. 

number 

estimate 

R 

R 2 

ID 

(2) 

(3) 

(4) 

1 

0.080 

0.939 

0.882 

2 

0.083 

0.945 

0.893 

FIG.  1. — Effects  of  Changes  in  Imperviousness  on  Characteristics  (TC  +  R)  and  (TC) 


If  development  is  predicted  within  one  of  the  drainage  basins  in  the  study 
area  and  if  drainage  area  and  slope  remain  relatively  constant  with  time,  then 
it  follows  from  Eqs.  I  and  2  that 


(TC  +  R), 
(TC  +  R)r 
(TC). 

and  - 

(TC), 


13) 

(4) 


hi  which  subscripts  e  and  f  refer  to  existing  and  future  conditions,  respectively. 
A  graph  of  the  left-hand  terms  of  Eqs.  3  and  4  versus  the  change  in  imperviousness. 
/,  is  shown  in  Fig.  1.  As  indicated  by  the  curves  and  indirectly  by  the  Tracor 
Report  (6).  Eqs.  I  and  2  are  significantly  limited  in  range,  particularly  for  practical 
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application.  Consider  as  first  approximations  the  following  two  examples. 

Example  I. — It  is  predicted  that  a  pristine  area.  (,  =  IY7r.  will  he  developed 
to  such  an  extent  that  at  some  time  in  the  future  it  will  he  considered  Hitt'  ; 
impervious.  In  cases  such  as  this  where  /,  -  I T4  initially  the  ratio  (/,//,) 
=  0  regardless  of  the  value  of  /(.  and  from  Fig.  I  or  Eq.  4.  <70,  0.  The 

preceding  results  are  impractical  even  for  very  small  values  of  (TC),. 

TABLE  2. — Results  of  Multiple  Regression  Analysis  in  Which  K  =  1.0  +  0.30  / 


-  - 

Coefficient 

Standard 

Correlation 

of 

Equation 

error  of 

coefficient. 

determination 

number 

estimate 

R 

R-’ 

(1) 

(2) 

(3) 

(4) 

8 

0.082 

0.937 

0.878 

9 

0.081 

0.948 

0.898 

FIG.  2. — Effects  of  Changes  in  Imperviousness  on  (TC  +  R)  and  (TC)  where  K  = 
1.0  +  0.30  f 


A  reasonable  estimate  of  (TC)f/(TC)r  for  a  condition  similar  to  the  preceding 
can  be  developed  from  the  following  formula  proposed  by  Kerby  (3): 


2  Ln 

3  V~S 


(5) 


in  which  t  =  the  time  of  concentration  for  overland  flow  within  a  catchment 
area,  in  minutes;  L  =  the  length  of  flow,  in  feet;  S  =  the  slope  of  the  surface, 
in  feet  per  foot;  and  n  =  a  retardance  coefficient.  In  a  situation  where  a  dense 
grass  covered  surface  will  be  completely  payed,  then 

. * 

(i),  \n,  / 


or 


0.18 


(7| 
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in  which  n,  and  »i,  are  0.02  and  0.80.  respectively. 

Example  2. — If  only  a  small  amount  of  development  is  predicted.  (K?  <  / 
<  5*/?.  for  a  relatively  pristine  drainage  area.  then,  in  most  cases,  it  is  reasonable 
to  assume  that  development  may  have  little  or  no  effect  on  the  time  of 
concentration.  Under  these  circumstances  (TO, /(TO,  should  he  equal  to  or 
nearly  equal  to  I . 

First  approximations  such  as  given  in  Examples  I  and  2  were  used  to  modify 
the  regression  expressions  (Eqs.  I  and  2)  as  originally  developed.  In  this  case 
various  transformations  were  tested  until  one  of  the  general  forms.  K  =  C, 
+  C,  /.  proved  applicable.  The  constants.  C ,  and  C\.  were  varied  until  the 
initial  approximations  were  reasonably  satisfied  and  an  optimum  degree  of 
correlation  obtained.  The  modified  relationships  (see  Table  2)  are 


(TC+  R)  =  17.01  K 

(TC)  =  11.54  K  0*' 


(8) 

(9) 


in  which  K  -  1.0  +  0.30  / 


(10) 


The  results  are  shown  in  Fig.  2  and  can  readily  be  compared  with  previous 
results.  The  ratios,  future  to  existing,  of  TC  +  R  and  TC  now  become 

( TC+R ),  /  1.0  +  0.30/  \ov> 

- L  =  ( - . (II) 

(TCfR),  V  1.0  +  0.30/,/ 

(TO,  /  1.0  +  0.30  I  y»*' 

and  - =( - )  . (12) 

(TO,  \  1.0 +  0.30/,/ 

From  Fig.  2  it  is  apparent  that  ( TC  +  R), /(TC  +  R)r  *  (TO, /(TO,  such 
that  a  practical  method  of  relating  TC  and  R  exists  for  a  particular  study  area. 
From  Eqs.  8  and  9 


TC 


TC  +  R 


11.54 


17.01 


/  DA  \": 


TC 


or 


TC  + 


/  DA  \n;4 

K‘i-r) 

/DA  Y,m 

—  =  0.68  K  00,  .  . (14) 

R  \  S  / 


For  practical  purposes.  Eq.  14  becomes 
TC 

- *0.68 . (15) 

TC+  R 

A  similar  analysis  yields 

R 

- *  0.32 . (16) 

TC+  R 


2 

\ 

* 

I 

< 

I 


I 

I 


1 

i 

I 


a 
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Relationships  presented  in  this  paper  can  he  used  ;is  a  guide  to  compute 
the  regional  unit  hydrograph  parameters.  TC'and  R.  for  existing  and  predicted 
values  of  imperviousness.  Modified  expressions,  such  as  those  developed,  are 
applicable  for  all  values  of  /  and  /,  /,. 
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